clear; clc; close all;

%% 固定参数
a0 = -2;
a1 = 0.83;
a2 = -0.34;
b0 = 1.3;
b1 = 0.1;
b2 = 1;
d1 = 0;
d2 = 0;

% 迭代设置
N_total = 500000;   % 总迭代步数
N_trans = 100;    % 舍弃瞬态步数
N_show  = N_total - N_trans;

% 初始条件
x0 = 0.1;
y0 = 0.1;
z0 = 0.1;

%% 分配存储数组
X = zeros(1, N_show);
Z = zeros(1, N_show);

%% 丢弃瞬态
x = x0; 
y = y0; 
z = z0;

% 先进行 N_trans 步迭代，不记录数据
for n = 1:N_trans
    [x, y, z] = vl_dmsin(a0, a1, a2, b0, b1, b2, d1, d2, x, y, z);
end

%% 记录后续 N_show 步
for n = 1:N_show
    [x, y, z] = vl_dmsin(a0, a1, a2, b0, b1, b2, d1, d2, x, y, z);
    X(n) = x;
    Z(n) = z;
end

%% 绘图
figure;
scatter(X, Z, 1, 'MarkerEdgeColor','none','MarkerFaceColor','b');
xlabel('x_n','FontSize',12);
ylabel('z_n','FontSize',12);
title('x, z相图','FontSize',14);
grid on;
